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Abstract 



Oscillations between ordinary and sterile neutrinos can generate large neu- 
trino asymmetries in the early universe. These asymmetries can significantly 
affect big bang nucleosynthesis (BBN) through modification of nuclear re- 
action rates. We study this phenomenon within a model consisting of the 
three ordinary neutrinos plus one sterile neutrino that can be motivated by 
the neutrino anomalies and the dark matter problem. We calculate how the 
lepton asymmetries produced evolve at temperatures where they impact on 
BBN. The effect of the asymmetries on primordial helium production is de- 
termined, leading to an effective number of neutrino flavours during BBN of 
either about 2.7 or 3.1 depending on the sign of the lepton asymmetry. 



I. INTRODUCTION 

Models with sterile neutrinos are currently favoured to accommodate all the neutrino 
oscillation measurements. It is difficult to simultaneously account for the solar, atmospheric 
and LSND anomalies with just three ordinary neutrinos and hence only two mass-squared 
differences. Therefore we require at least one additional neutrino if we wish to explain all 
three of the observed neutrino anomalies. Constraints from the Z-width suggest that any 
additional neutrinos must be sterile. Oscillations between ordinary and sterile neutrinos in 
the early universe have been studied in and found to generate large neutrino asymme- 
tries ||||, which can significantly affect big bang nucleosynthesis (BBN) reaction rates, and 
can result in an effective neutrino number of less than three during BBN || . In this paper, 
the techniques developed in Ref. M will be applied to the four neutrino model of Ref. M 
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which has almost degenerate and v T . This model is distinct from the four neutrino model 
of Ref. || where a mass hierarchy between and v T was assumed. 

A motivation for studying how neutrino oscillations can affect N***, the effective neutrino 
number during nucleosynthesis, arises from discrepancies between the observed deuterium 
abundance and standard BBN. The primordial D/H ratio can be used to give a sensitive 
determination of the baryon to photon ratio rj which, given the inferred primordial 4 He mass 
fraction, can be used to predict the number of light neutrinos N%*' . 

Although there are numerous measurements of the deuterium abundance, obtaining an 
upper bound to the abundance is quite difficult, and there are conflicting observations. A 
high deuterium measurement of D/H = (1.9 ± 0.4) x 10 -4 was found in Ref. ||, leading to 
r] ~ 2 x KT 10 , while in Ref. § a low deuterium result of D/H = (2.3 ± 0.3 ± 0.3) x 10~ 5 
was obtained, suggesting that r\ ~ 7 x 10 -10 . From these values of rj, and the helium mass 
fraction Yp, the effective number of neutrinos can be calculated. By way of example, the 
results 

N e J f = 2.9 ±0.3 high D/H, 

N e J f = 1.9 ± 0.3 low D/H, (1) 



were obtained in Ref. |ItJ. (For other calculations of N„** see Ref. |n|.) Hence, if the low 
deuterium result is correct, standard BBN would require an effective neutrino number of 
less than three (provided that the value of Yp used in Ref. jnj is correct). Clearly N*" < 3 
would imply new physics. Further observations are required to pin down the D/H value 
(and also Y P ) although it now seems as though the low D/H result is favoured P,|T2 . 



One possibility for the new physics is to look at the neutrino sector and study mechanisms 
by which the effective number of neutrinos can be made less than three. A possible way 
of doing this is to create an electron-neutrino asymmetry, as this would directly affect the 
reaction rates which determine the n/p ratio just before nucleosynthesis. The n/p ratio 
controls the helium mass fraction Yp via the equation 

= _ A(n _> p)Yp + A(p - n)(2 - Y P ), (2) 

where the reaction rates 

A(n — * p) ~ X(n + v e — > p + e~) + A(n + e + — > p + P e ), 

\(p — > n) ~ A(p + e _ — > n + v e ) + \{p + v e — > n + e + ), (3) 

depend on the momentum distributions of the species involved. The processes in Eq.([3]) 
for determining n p are valid for temperatures above about 0.4 MeV, below which the 
weak interaction rates become frozen out and neutron decay becomes the dominant factor 
affecting the n/p ratio. An excess of v e over v e would reduce the n/p ratio, due to the 
modification of the neutrino momentum distributions through non-zero chemical potentials, 
thereby changing the rates for the processes in Eq.(|3]). Neutron decay is not significantly 
altered by lepton asymmetries. A change in the helium mass fraction Yp can be related to 
a change in the effective number of neutrinos by [13| 



5Y P ~ 0.012 x 5N e J f . (4) 



2 



A large neutrino asymmetry could simply be postulated, but it is more interesting to look 
at physics of the neutrino sector which will allow a large neutrino asymmetry to evolve as a 
dynamical variable. 

A four- neutrino model was studied in Ref . , involving the three ordinary neutrinos and 
one sterile neutrino, where the mass hierarchy m Ur ^> m u „, rn Ue , m Vs was considered by way 
of example. In this case oscillations between u T and u s resulted in an excess of v T over v T 
thereby generating a large tau-neutrino asymmetry and hence a large tau-lepton number. 
v T — u e oscillations also occurred, and it was shown that this allowed some of the tau-neutrino 
asymmetry to be transferred to an electron-neutrino asymmetry. The effective number of 
neutrinos found in Ref. || was either 2.5 or 3.4. The two values arise because there is an 
ambiguity involving the sign of the asymmetry and hence the prediction for (see Ref. 

for a discussion of this issue). For a positive asymmetry, 5N*ff ~ —0.5 was obtained over 
a range of mass differences 5m^ s ~ 10 — 3000eV 2 , while for a negative asymmetry the result 
was 8N*ff ~ +0.4. So we have the interesting situation that N^f < 3 can be achieved 
although the precise value of N*'* is model dependent as this paper will illustrate. If future 
measurements can pin down N%** precisely enough, then this quantity can be used to help 
discriminate various competing models explaining the observed neutrino physics anomalies. 



II. MODEL OF NEUTRINO MASSES 



The model we will consider in this paper is again a system of four neutrinos: u T , u^, u e 
and u s . In this model the v e and v s are assumed to be very light (<^C leV), and v T and are 
taken to have nearly degenerate masses. This model has been proposed in M as a pattern 
of masses that will allow the data from the solar, atmospheric, and LSND experiments to 
be satisfied as well as being consistent with dark matter models. 

In this scenario, the solar neutrino problem is explained by the small angle MSW solution 

" 5 eV 2 and sin 2 26 
oscillations 



with oscillations between v P and z/, with <5m 2 D ~ 10" 



10~ 2 14 



atmospheric neutrino deficit is explained by v T 



10 3 eV 2 and sin 2 26 ~ 1. The LSND data requires a 
0.2 - 10eV 2 and sin 2 26 ~ 3 x 10~ 2 - 10~ 3 



TBI with (5m 2 , ~ 



while the 

io- 2 - 



' e solution with <5m 2 e in the range 
T6f . It has been argued [0 from considerations 



of structure formation in the early Universe that the optimal value for the degenerate 
and v T masses is 2.4eV. Note, however, that the creation of a large L M and L T neutrino 
asymmetry should impact on the favoured neutrino mass for structure formation and dark 
matter. For this reason it may well be necessary for the optimal mass of approximately 2.4 
eV to be recalculated, taking into account the modified neutrino number densities. 



III. NEUTRINO ASYMMETRIES AND CONSEQUENCES FOR BBN 

For oscillations between v a and u s , the weak eigenstates are linear combinations of the 
mass eigenstates V\ and u 2 , 

v a = cos 6*o z/i + sin #0^2, 
u a = — sin 9qV\ + cos #0^2, (5) 
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where 8q is the vacuum mixing angle (we assume that 9q is defined such that cos2#o > 0). 



The mass-squared difference between the two eigenstates is defined as 8m c 



rat, 



■mf. The 



neutrinos interact with matter, so a matter mixing angle 9 m is defined, which is related to 
#o by 



sin 2 20„ 



sin 2 26n 



sin 2 26> + (b ± a - cos 26» c 



^2' 



(6) 



The term (6±a) is related to the effective potential due to interactions of the neutrinos with 
matter: 



V 



2p 



±a), 



(7) 



where -/+ corresponds to vjv and p is the neutrino momentum. The functions a and b are 
given by 



4:V2((3)G F T 3 L^p 



ir 2 5m 2 „ c 



4V2C(3)C F r 4 A 



p 



ir 2 8m 2 as M^ 



(8) 



where A e ~ 17 and A T ^ ~ 4.9 (TP 



The function is defined as 



L (a) = L a + L T + + L e + 77, 



(9) 



where r\ is approximately equal to the baryon to photon ratio, and the L u 's are the lepton 
asymmetries which are given by 



n„„ - n,- 



(10) 



where n v and n 7 are the neutrino and photon number densities, respectively. 

The neutrinos have a distribution of momenta and only neutrinos with momentum of a 
certain value will be resonant at a given time. The momentum of the resonance will move 
through the neutrino distribution as the temperature decreases, creating a neutrino asym- 
metry at the particular value of the resonance momentum at that temperature. The evolving 
neutrino asymmetry also non-linearly affects the evolution of the resonance momentum. The 
asymmetry is distributed across the momentum distribution as weak interactions keep the 
distributions in thermal equilibrium. 

The resonance occurs for 



(b±a) = cos2# ^ 1, (11) 

for sin 2 2^o <^ 1- The a and b terms dominate this condition at different temperatures. We 
may obtain a qualitative understanding of the behaviour by neglecting the thermal spread 
of momentum and using the average value (p) ~ 3.15T. At high temperatures the (b) 
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term dominates, and in [|J it was found that for (b) < 1 oscillations create lepton number 
while for (b) > 1 the oscillations are lepton number destroying. The dependence of (b) on 
temperature is (b) ~ T 6 , so when the temperature drops to where (b) ~ 1, which we denote 
by Tc, the oscillations creating lepton asymmetry dominate and lead to an exponential 
growth in lepton number. In Ref. Q it was shown that for 5m 2 < (and \5m 2 \ ^ 10~ 4 eV 2 ) 
large lepton asymmetries are created provided that 

/ e v 2 \ 1/6 

*'» * 5 x 10- (^J . (12) 

This result is independent of the initial value of the asymmetry, provided that \L a \ ^ 10~ 5 
Q. For 5m 2 in the range 1 — 10eV 2 , we have Tc — 16 — 25MeV, obtained from the condition 
(b) = 1. The rapid growth drops off as the temperature decreases and when the temperature 
has decreased to below about Tq/2, the (a) term has become larger and we can now neglect 
the (b) term so that the resonance condition is approximately given by 



(a) ~ cos 20 o ^ 1- (13) 

Collisions and oscillations both affect the generation of lepton asymmetries. At high tem- 
peratures (T > Tc/2) collisions affect neutrinos to a greater extent than oscillations, and 
these in turn affect lepton number because the z/'s and z/'s interact at different rates in a 
CP asymmetric background medium. (A CP asymmetric medium is required by the cosmo- 
logical baryon-antibaryon asymmetry.) At lower temperatures, the collision rates are not 
as rapid and oscillations become coherent, so that oscillations take over from collisions to 
control the growth of lepton number. Thus, for T < Tq/2, MSW |18[ transitions become 
important. 

We have two oscillation resonances to consider, z/ TjM — z/ s , which will generate a large 
Lt./j, asymmetry and z/ T)M — z/ e , which will then transform some of the asymmetry into the v e 
sector to produce L e . Oscillations between v T and z/ M or z/ e and v s may be neglected as they 
have much smaller mass-squared differences. 

The asymmetries are calculated, starting from an initial temperature of T = Tq/2 where 
\L T \ <C 1 and T is low enough for oscillations and MSW transitions to dominate over 
collisions. The initial number density of sterile neutrinos is taken to be negligible (which 
will occur for a range of parameters) in accord with the results of the numerical calculation in 
Ref. 0] for the initial creation of lepton number at T = Tc- We study how the asymmetries 
evolve down to temperatures where they become important for BBN. Note that for the 
moment we will assume that L T > for definiteness. 

At T > Tc where the b term dominates, v T —>■ v s and v T — > v s resonances both occur at 
the same momentum, but as the temperature falls to T ~ Tc/2 the v T — > v s and v T — > v s 
oscillations no longer have the same resonance momenta. This allows the asymmetry to 
grow very large because when most of the z/'s have approximately the right momentum to 
be near a resonance, the bulk of the z/'s are far away from a resonance. It was found in 
that if the created L a is positive the antineutrino resonance momentum has generally low 
values Pres/T ~0.2 — 0.8atT~Tc</2 while the neutrino resonance momentum is toward the 
tail of the neutrino momentum distribution (see Fig. 2 of Ref M for an illustration). Hence 
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the neutrino resonance is neglected and we will only be concerned with the antineutrino 
resonance. 

We will neglect the v T — mass difference and thus make the approximation that the 
v T — v s and — v s oscillations have exactly the same resonance momentum, pi, and the 
v T ,ii — v e resonances have momentum p2- These are obtained from the resonance condition 

a ~ cos 2^o — 1, (14) 
where a is given by eq.(§). The resonance momenta are given by 

ir 2 5m 2 1 



Pi 



K 2 $m 2 1 n _v 

P2 ~ 4v^C(3)G F T3(Z>) 1 J 

Since v T and are nearly degenerate and their matter effects are identical, the calcula- 
tion is symmetric with respect to and v T) so L T = L^. We therefore find that 



L (t) = 3L T + L e ~ 3L 7 
L (t) - L (e) = L T - L e ~L T . (16) 

This means that the v T ^ — v e resonance occurs for a momentum which is approximately a 
factor of three higher than for the v T M — v s resonance, 

P2 ^ 3pi- (17) 

Because the v T — v s and — v s oscillations have the same resonance momenta, instead 
of having two separate resonances, the v T — v s and — v s resonances completely overlap. 
Therefore we must consider the three-flavour system consisting of z/ r , and v s . To take 
into account that there is only one sterile neutrino state for both v T and z/ M to couple to, 
we may form two linear combinations of v T and v^. One of these oscillates with i> s , and 
undergoes MSW conversion at the resonance, while the other v T j linear combination does 
not oscillate with v s . This state is decoupled, and is not converted into the sterile state at 
the resonance. Our equations include factors of \ to account for this situation. The mass 
eigenstate which is dominantly a v s at high density, becomes a Vrjv^ linear combination 
at low density, so that a v s which goes through the resonance has a 50% chance of being 
converted to v T and a 50% chance of being converted to v^. Similarly, since one linear 
combination of v T and only will be converted to a v s upon passing through the resonance, 
both v T or have only a 50% probability of being converted to the sterile state. 

Naively one might believe that by making and v T degenerate, so as to have two heavier 
neutrinos oscillating with v s instead of only one heavier neutrino oscillating with the sterile 
neutrino, that twice the electron-neutrino asymmetry would be produced. However, this 
will not occur, because there is only one sterile eigenstate for the two linear combinations 
of v T and to oscillate with. Due to the factors of \ arising from this three-flavour nature, 
we would not expect the electron-neutrino asymmetry created to be any larger than for the 
model of Ref . || , where only a heavy v T oscillates with v s (in fact it turns out to be smaller 
-see later). 
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Numerically integrating the quantum kinetic equations, we find that the adiabatic limit 
of the MSW effect holds provided 



sin 2 26 n > 10 



-9 



10 



-10 



(18) 



for Sm 2 ~ 10eV 2 . The adiabatic limit implies all neutrinos that pass through the resonance 
will undergo MSW conversion. This means the lepton number created is basically equal to 
the number of antineutrinos minus the number of sterile antineutrinos that pass through 
the resonance, and the growth of lepton number is governed by how quickly the resonance 
momentum moves through the neutrino momentum distribution. 
The rate of change of the lepton numbers is therefore given by 



dL,. 



dT 



d 


(Pi 


df 


\T 




dL Uc 




dT 



Xo 



2X, 



d 




dT 


m 


d 




df 





(19) 
(20) 



where Xi(X2) expresses the difference between the number of v T and u s (u e ) with the right 
momentum to pass through the resonance and is given by 



X x 



X, 



1 T 

2 ra 7 

1 T 

2 7i.v 



(N, T - N Ps ) 



(N Pt - N D 



p=pi 



P=P2 



The factor of 2 in Eq.(pO|) and the factors of | in Eqns. 
nature of the problem as discussed above. 



(21) 



arise from the three-flavour 



The N's are the momentum distributions of the neutrinos, so that 



n = J Ndp, 

which in thermal equilibrium are just given by a Fermi distribution 



(22) 



N 9 



2tt 2 i 



exp 



P+P-v 

T 



(23) 



_d_ 

dT 



( Pres \ 
t J 



terms describe how quickly 
because we 



where fip is the antineutrino chemical potential. The T 

the resonance momenta move, and are expressed in terms of the variable r 
want to distinguish the roles of expansion and oscillations, and unlike p, the quantity ^ 
does not change due to the direct effect of the expansion of the universe. Evaluating the 
derivatives, we obtain 



dL, 



dT 
dL,, 



2D A + CE 
A + E(A — B) ' 
C B dU 



dT 



A + A dT 



(24) 
(25) 
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where 



TL( T ) J ' 2 TLW ' 

_ AX lPl AX 2 p 2 

E " T LW-IW (26) 

In deriving eqns.( |2^ , |23| ) we have used that ^ (^) < ^, because the resonances (which 
start at low values of approximately pi/T ~ 0.2 — 0.8 and p 2 /T ~ 3px/T) move through the 
momentum distributions to higher values as the temperature decreases. 

Numerically integrating eqns. (p4| , 25) , we find that the final lepton asymmetries generated 



are 

Ljh = L T /h ~ 0.16, L e /h ~ 6.7 x 10" 3 , (27) 



for 5m 2 in the range 6-8eV 2 , where h = /T 3 . 

The L e asymmetry then affects BBN rates via the modification of the neutrino momen- 
tum distributions. Weak interactions keep the neutrino distribution in thermal equilibrium 
for T > lMeV, however for 5m 2 ~ 10eV 2 the neutrino asymmetries only become large at 
temperatures where the the distributions are beginning to go out of equilibrium. Because 
the neutrino distributions are not completely thermalised, we cannot simply describe them 
in terms of chemical potentials, so for this reason we keep track of the neutrino number 
densities numerically in momentum cells. We assume complete MSW conversion at the res- 
onance, which is a good approximation for the large range of parameters given by Eq. (|TS|) . 
so that numbers densities of sterile and active neutrinos are interchanged at the resonance 
momentum such that 

n*M - + tfpJCPi) = 

N*,,M - \{N>M + N, T Jpi)), (28) 
and similarly for the v T[i — v e resonance. Cells are refilled according to the interaction rates 

d ( -^actual \ , ,/ \ ( ^ equilib ^actual \ / nn \ 

It [—) = r(p) {— ivTV • (29) 

where 

1 v 2 

N = — v . . . (30) 

and N equi i ib is given by the equilibrium distribution Eq . (p3|) . The chemical potential is 
obtained from the lepton number via the equation 



L ~ 



24C(3) 



^{fi v - }i p ) - 6(fii - /2?) In 2 + {fit - /i 3 ) (31) 
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where jli = m/T. Equation (|31~D is an exact equation for p, v = —flu, otherwise it holds to 
a good approximation provided that jli ^ 1. For 5m 2 ~ 8 eV 2 significant lepton number 
(here significant means larger than about 0.01) is not created until T < 2 — 3MeV. At these 
temperatures, the neutrinos and anti-neutrinos have already chemically decoupled. (Note 
that Tjg C ~ 3 MeV and T^ T C ~ 5 MeV are the chemical decoupling temperatures). Because 
of this, \x v ~ 0, while the anti-neutrino chemical potential ftp continues increasing as per 
eq.©. 

The interaction rate T(p) in eq.(|29|) is 

r W - ^> (32) 
where (r) is the thermal average of the total neutrino collision rate given by HJ 

(T Ue ) = 4.0G 2 T 5 , 

(F„ T J = 2.9G 2 F T 5 . (33) 

Integrating the rate equations for the processes given in eq.(H), using the modified neu- 
trino momentum distributions (N actua i(p)) , the change in the neutron to proton ratio and 
hence Yp is determined. We find the size of the effect in this case is given by 5Yp ~ —0.0023 
which corresponds to a change in the effective number of neutrinos of 5N^f ~ —0.19. 

Excitation of the sterile neutrino, and the modified momentum distribution, will change 
the energy density of the universe, and this will also change the effective number of neutrinos. 
The number density and average energy of the neutrinos were calculated from the modified 
momentum distributions. For 5m 2 in the range 6 — 8eV 2 we get 

— = 0.99, Vt^ = 0.99, 

n Q (p) 

n p rip (E p ) (Ep ) 

^ = ^ = 0.53, Vt^ = Vt^ = 0-91, 

n n (p) (p) 

1^ = 0.91, ^1 = 0.97, (34) 
n (p) 

where n is the number density with zero chemical potential and (p) = 3.15T is the average 
momentum for a Fermi-Dirac distribution. (The chemical potentials for the neutrinos are 
very small and so the distributions are not significantly modified.) These changes in the 
energy density correspond to 8N%" ~ —0.09. As we have an ambiguity concerning the 
sign of the lepton asymmetries and hence also the sign of 8N„*f due to modification of the 
nuclear reaction rates, we have two possible values for the overall change in the effective 
number of neutrinos. 

In fig.(|l]) 5N*ff is plotted as a function of the mass squared difference. For the mass 
range 5m 2 ~ 6 — 8eV 2 , we find: 

8NI** ~ —0.3, positive asymmetry, 

8N„* ~ +0.1, negative asymmetry, (35) 
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This remains approximately constant to above where the NOMAD experiment cuts off the 
LSND data at Sm 2 ~ 10eV 2 Eofl, while for small Sm 2 we have 5N?/ f -> 0. 



The L e asymmetry and hence the 5N^f obtained is smaller than that obtained for the 
model of Ref . . The L T ,fi asymmetry created in the present model is roughly half as big as 
the L T asymmetry obtained in Ref. ||, because only one of the the two linear combinations 
of U/j, and v T oscillates with u s , as discussed above. Naively one might expect that with both 
v T — v e and v^ — v e oscillations producing L e that an effect of comparable size to Ref. || would 
be obtained, however, this does not turn out to be correct. The crucial difference between 
the two cases is the separation of the two resonances. The separation of the resonances is 
given by Eq. ( |T7| ) , and whereas we have p<i — 3pi , in the previous case the resonances were 
closer together with p2 — 2p±. This is important because the resonance which creates L e 
is at a higher momentum than the resonance which creates L T , so that by the time L T has 
grown sufficiently large to start producing significant L e , the p^ resonance has already moved 
a considerable way through the neutrino momentum distribution. For the present model, 
P2/T has already moved past the average value of p/T ~ 3.15 and is out toward the tail of 
the distribution, where the number density is much lower, so that the L T asymmetry cannot 
be transferred to L e as efficiently. 



IV. CONCLUSION 

We conclude that the neutrino oscillations in the four neutrino model of Ref. [0] can 
modify the effective number of neutrinos during nucleosynthesis and give Nj;'* < 3. How- 
ever, it turns out that the effect is relatively small, SN^'* ~ —0.3 (if L e > 0). If future 
measurements confirm the low D/H result as well as the Yp value used in Ref. [[10[| , then this 
model will become cosmologically disfavoured. On the other hand, if future measurements 
favour N*" ~ 3 then this model will remain viable. 

There is also the interesting possibility that the neutrino mass favoured for dark matter 
fT/| will have to be recalculated. In fact, from Eq.(^) the number of heavy neutrinos is 
actually about 1.5 because about 50% of the u T and have been converted into light v s . 
This would suggest that the optimal mass for hot+cold dark matter senarios is somewhat 
heavier, about m UT = m u = 3.1eV rather than 2.4eV, for this particular model. This implies 
that the 5m 2 relavant for LSND is about 9eV 2 rather than 6eV 2 [ pTf . 
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FIGURES 

FIG. 1. Change in the effective number of neutrinos as a function of the mass-squared difference. 
The solid line corresponds to a positive electron neutrino asymmetry, and the dashed line a negative 
asymmetry. 
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